Significance of the hyperfine interactions in the phase diagram of LiHo x Yi_ x F 4 
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We consider the quantum magnet LiHo x Yi_ x F4 at x = 0.167. Experimentally the spin glass to 
paramagnet transition in this system was studied as a function of the transverse magnetic field and 
temperature, showing peculiar features: for example (i) the spin glass order is destroyed much faster 
by thermal fluctuations than by the transverse field; and (ii) the cusp in the nonlinear susceptibility 
signaling the glass state decreases in size at lower temperature. Here we show that the hyperfine 
interactions of the Ho atom must dominate in this system, and that along with the transverse 
inter-Ho dipolar interactions they dictate the structure of the phase diagram. The experimental 
observations are shown to be natural consequences of this. 



LiHo x Y!_ x F4 is widely considered to be a model quan- 
tum Ising magnet. The strong easy axis crystal field gives 
the large spin Ho ion an Ising character, behaving as a 
2-level system when temperature T < 10 K; quantum 
fluctuations between the 2 Ising states are induced by a 
tunable transverse magnetic field H±. If one then ne- 
glects both hyperfine (hf) and transverse dipolar (TD) 
interactions, the system is described by a transverse field 
Ising model (TFIM). Experimentally, for x = 1 the sys- 
tem orders ferromagnetically at low T and when 
x < 1 the longitudinal dipolar interactions become ran- 
dom and frustrated, and for x = 0.167 it was observed 
that the system orders in a spin glass (SG) statejj 0] • 
Much of the recent theoretical interest in quantum spin 
glasses 0,0,0 has been driven by these experiments^,]]!, 
and both general theories on quantum spin glasses (see 
e.g. Ref. 0) and specific models of the quantum Ising 
SGH,|i| are commonly checked against their results. 

However, for x — 0.167 the hf interaction is larger than 
the mean dipolar interaction, and in part (i) of the Let- 
ter we show that it actually dominates the low-energy 
physics in this system. The longitudinal hf term forces 
the two ground doublet states to have a definite, and 
opposite, nuclear spin|lOl 111) . Low transverse magnetic 
fields (hbH± <C f2 j where fl is the gap to the first 
excited electronic state) cannot then couple or induce 
quantum fluctuations between these Ising- like states, but 
only renormalize their effective spins (and hence the ef- 
fective dipolar interaction). Transverse hf interactions 
can induce fluctuations between these states; however, 
when [abH± <C f2 Q we see that these fluctuations are 
very weak. Therefore, at low transverse fields, the low- 
energy effective Hamiltonian is given by the classical Ising 
model, with renormalized parameters. We thus explain 
why gusHc ^> kT c in the experiment i.e., why the 
transverse field coupling required to destroy the glass is 
so much larger than the thermal energy required to do 
the same: While T c at zero field is given by the mean 
dipolar interaction, the temperature dependence of the 
transition line on H± is dictated by the larger hf inter- 
action that controls the renormalization of the effective 
spin, and the critical transverse field is controlled by the 
even larger energy scale of ft . 



To completely understand how fluctuations are 
switched on by transverse fields, we show that TD terms 
must be included. By doing this in a simple way we ob- 
tain a satisfactory quantitative agreement with the ex- 
perimental phase diagram. We then consider the dimin- 
ishing of the cusp in the nonlinear susceptibility ^3, found 
in the experiment |3j as temperature is reduced. A natu- 
ral consequence of our theory is that this may be a result 
of the renormalization of the effective spin, rather than 
a first order phase transition. 

The Hamiltonian describing the LiHo x Yi_ x F4 system 
is given by a sum of crystal field [ToLIT^ . Zeeman, hf, and 
inter-Ho interaction terms: 

H = H c{ + H Z + H M + H mt . (1) 

The Ho atom has a total angular momentum of J = 8, 
and nuclear spin I = 7/2. The Zeeman term Hz = 
~ Si QjI^bH ■ J h and H hi = Aj J2i h ' Ji is tnc hf in- 
teraction. Hi n t — U°f ' JiJj is dominated by the 
dipolar interaction|12|. The interaction of the external 
magnetic field with the nuclear spins is small, and is ne- 
glected here. 

The crystal field H c f leaves only three electronic states 
at low energy; a low-energy doublet, denoted |f) and |J.) 
here, and an excited state |r^) roughly fi G = 10.5 K 
above these. Because of the strong crystal anisotropy the 
TD and hf terms are usually neglected. Neglecting the 
longitudinal hf interactions as well, and in the presence 
of transverse magnetic field H± , the low-energy effective 
Hamiltonian for the system would then be the simple 
TFIM0: 

H = -T, V ^ T ^ - A °™ E ^ - w 

where fj is a Pauli vector describing the two-level ef- 
fective electronic spin at spatial position r = Yj , the j-th 
site, V** w 150J7?? (with |V-f | ~ 0.3 K for nearest- 
neighbour Ho spins), and A Q is the transverse field- 
induced coupling between the eigenstates ||) and |J.) 
of t z . For small fields A D oc H\ coming from the 
second-order coupling, mainly via the state IT^) [so that 
A Q ~ 9([1bH±) 2 /0 o ]. In this simple model quantum 
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FIG. 1: Splitting of the electronic low-energy doublet ( j and 
J.) due to the longitudinal hyperfine interaction. The ground 
state doublet, a and a have a definite and opposite nuclear 
spin, ±7/2. Transverse magnetic field couples states with the 
same nuclear spin, as is shown by the dashed lines. 



fluctuations are already important at fairly small trans- 
verse fields, and we expect a T = quantum phase tran- 
sition (QPT) when A - V a , where V D = (£\ V t f) is the 
average dipolar interaction. 

For the undiluted system, x = 1, it was indeed 
shownQ, IT2I fl3| that the ferromagnetic paramagnetic 
(PM) transition line is well described by this model, with 
a low temperature correction due to hf interactions. How- 
ever, in the dilute case V D is reduced by a factor x, so the 
hf interaction is much larger than V . The model is 
then inadequate, and we must include the hf interactions 
from the beginning. We now introduce the interactions 
neglected in the model J5J) consecutively, emphasizing the 
significance of each on the structure of the phase diagram. 

(i) Hyperfine Interactions: We first consider the regime 
(j.bH± <C ri in which the transverse hf interactions are 
negligible. The longitudinal hf interaction H£ { — AjI z J z 
splits each of the states |f), |J.) into an eightfold multiplet 
of nearly equidistant levels, with separation ~ 205mK[l(J 
between adjacent levels (Fig|TJ. The two lowest energy 
Ising states in zero field are now a = If, —7/2) and a = 
||, 7/2), and when H±_ 7^ 0, these become 
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with coefficients c\ = aA , and C2 = a[A 
where a = [A 2 , + (A - A 2 + A 2 ,) 2 ]- 1 / 2 and 2A « 1.4 
K is the energy difference between the states a(a) and 
6=|i,-7/2)(6=|T ) 7/2)). 

The crucial point here is that the transverse magnetic 
field H± does not couple the relevant Ising doublet states: 
a and a (it actually couples a to b and a to b). This im- 
mediately invalidates the TFIM J2Jl, since H± no longer 
induces fluctuations. The only effect of the transverse 
field is to re-normalize the spin: For the state | +) one 
finds 

«)=r / «(0)) ; r^( C 2 -c 2 ), (4) 

and (erf.) = —(erf). One can then absorb the renor- 
malization into the dipolar interaction, and therefore our 
system reduces at low transverse fields to the classical 
Ising model: 
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(5) 



where s is a spin half matrix in the space of the states 
I +) and I -) and Vg* = rfVfi* . 

As a result, the SG paramagnet transition line is given 
(in this model, neglecting the later to be shown signifi- 
cant TD interactions) by the relation T c = n 2 T c (0), and 
the dependence on the transverse field is through i](A Q ). 
Importantly, this result is independent of the approxima- 
tion used to treat the longitudinal dipolar interaction, 
and relies only on the assertion that T c scales linearly 
with the interaction. For A D <C A one has from Q that 

A 2 

t] = 1 — an( l defining e = (T c — T)/T c one finds that 
for A /A, e < 1 (i.e., small H± and T — T c ) the phase 
transition line A C (T) obeys the relation A c = A^fe. For 
the initial TFIM one would instead get A c ~ V Q y/e. 
The difference arises because the hf energy scale A dic- 
tates the reduction of the effective spin (and equivalently 
the effective dipolar interaction). At T — T c /2 one finds 
that A c w A, and for A a /A > 1 we have ?/ = A/A a , so 
that A c = Ay/ (V /T). Thus, as long as the transverse 
hyperfine interactions are negligible ({1bH± fl ) there 
is no QPT. 

We now turn to the discussion of quantum fluctuations 
in this system, and we thus consider the transverse hf 
coupling = Aj(I+J- + J_ J+)/2. This term changes 
the z component of the nuclear spin by coupling the elec- 
tronic states |t) and ||) with (mainly) the state | T l 2 ) at 
energy f2 D , thus inducing fluctuations between | +) and 
I — ). However, unless £ = (||/Xb^± | r|) ~ (cor- 
responding to a field H± ~ 2T), these fluctuations are 
much smaller than V , and cannot induce a QPT. To 
see this more quantitatively we diagonalize the single Ho 
Hamiltonian in a transverse magnetic field, i.e., we diag- 
onalize H = H c f + Hz + Hm, in the 136 eigenfunction 
space (17 crystal field states * 8 nuclear states) using the 
parameters used in Ref. [l^. We then plot in Fig. 2 the 
results for the splitting A(H±) between the two lowest 
levels induced by the combination of H± and H^. For 
H± < 3T the two lowest levels are well separated from 
the higher energy levels, and we can replace the classical 
Ising Hamiltonian (J5J by an effective Hamiltonian 



A(H x ,A,n o )J2s X i 



(6) 



We then expect a zero temperature QPT when A(H±) — > 
A C (H±), such that A c (Hj_) m V (Hx). From Fig. 2 we 
confirm that this happens when H± ~ 2T. Note there are 
now three significant energy scales in the problem. T c is 
dictated by V , while A c is related to both A, which dic- 
tates the renormalization of the effective dipolar interac- 
tion, and to f2 . Interestingly, for any practical dilution x 
the lower bound for the critical transverse magnetic field 
is set by f2 , since only at £ « f2 Q do quantum fluctuations 
between the relevant Ising states become appreciable. 

We may summarize our analysis so far by saying 
that the longitudinal hf interactions invalidate the usual 
quantum Ising model for the LiHo x Yi_ x F4 system, in- 
stead creating a classical Ising model whose energetics 
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is determined by the longitudinal dipolar interactions, 
renormalized by longitudinal hf interactions. Adding in 
the transverse hf interactions brings back an effective 
transverse field term A(H±) into the effective Hamilto- 
nian, but A(H±) does not become important until much 
higher transverse fields than in the original model (and 
it switches on with field in a quite different way from 
the usual form A Q oc Hj_). We emphasize that A should 
be used instead of A D in analysis of the phase transition 
and, in particular, for the determination of the critical 
exponents in the system. 

Can our simple model approximation explain the ex- 
perimental phase diagram? Since according to our anal- 
ysis up to now, the position of the phase transition line 
is mainly dictated by single atom properties and the hf 
interaction, we now derive a phase diagram treating the 
single Ho exactly, and the inter-Ho interactions by us- 
ing a mean-field (MF) approximation in which each spin 
feels an interaction strength proportional to its average 
J z , i.e., using the MF Hamiltonian 

Hmf = H c [ — g,j^BH±Jf 

3 

+ .1,^ /;../; ■ Y^\„ ./; (7) 
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In Fig. |3| we plot by dashed line the SG-PM transition 
obtained within this approximation. V is fixed by the 
experimental value of T c at zero transverse field. In or- 
der to compare our results to Fig. 1 of Ref. Q we plot 
the phase diagram as function of T and A D , which is 
given by the splitting of the electronic levels when the hf 
interactions are neglected. 

This mean field result does capture two main features 
of the experiments, viz. (i) the much larger energy associ- 
ated with quantum disordering of the SG state than with 
temperature disordering, and (ii) the roughly straight 
transition line in the A Q ,T plane in most of the param- 
eter regime. However two noticeable differences are ap- 
parent, viz., (i) the experimental T = quantum critical 
point occurs at lower transverse field and (ii) the shape of 
the experimental transition line near the zero field tran- 
sition point is quite different from the predicted square 
root behavior. These differences are not a result of our 
MF approximation. As discussed above, the square root 
behavior at low transverse fields is a direct consequence 
of the renormalization of the spin by the transverse mag- 
netic field, and is a single atom property. Furthermore, 
the experimental graph suggests that the quantum fluc- 
tuations induced by the transverse hf interactions are 
significant already at fields of the order of 1 Tesla, in 
contrast with FigEl 

(ii) Transverse Dipolar Interactions: The differences 
just noted suggest that the TD interactions [e.g., (a/?) = 
(zx)), neglected in the Hamiltonian J7J), are also signif- 
icant in the diluted LiHo x Yi_ x F4 system. Unlike for 
x = 1, where the transverse terms cancel by symme- 
try, in the diluted system this is not the case- the TD 
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FIG. 2: The splitting A(H±, A, fl ) between the ground and 
first excited state of the single-ion Hamiltonian as a function 
of H± (solid line) calculated by exact diagonalization. This 
measures the quantum fluctuations between the states | +) 
and | — ). The dashed line is the averaged longitudinal dipo- 
lar interaction, taking V = 0.13 for H± = to match the 
experimental T c at zero field. 

terms add an effective magnetic field at each Ho loca- 
tion. This field adds both an effective random longitu- 
dinal term at each site|l4|. and a transverse term which 
can induce fluctuations even in the absence of H±. Wc 
consider the latter effect here, and the former in con- 
nection with the nature of the phase transition below. 
For [1bH± « V or larger it would seem that the in- 
duced transverse fields are random and could enhance or 
decrease the effective transverse field at each location. 
However, the external magnetic field breaks time rever- 
sal symmetry, and the configurations in which the TD 
interactions enhance the effective magnetic field at the 
Ho sites are energetically favorable, as can easily be seen 
by comparing the two Ho states ||, — >) and \l, —*). Note, 
that local strain fields result in random effective fields in 
the transverse direction as well. However, their magni- 
tude can be estimated to be smaller than that of the TD 
interactions ^3 j an d more importantly, they are symmet- 
rically distributed with zero mean. Their affect on the 
effective transverse field is therefore neglected here. 

In order to demonstrate the significance of the TD in- 
teraction, we now redo the MF calculation with the sim- 
plified approximation of taking the effective transverse 
field as H± = H± + Hf . The effective field due to the 
dipolar interaction Hf is the only free parameter in our 
calculation. The solid line in Fig. [21 is a plot of the SG- 
PM phase transition line including the TD interaction, 
with Hf = 0.65T.[15j The experimental transition line is 
reproduced satisfactorily, including its shape near T c and 
the critical transverse field at low T, using a value of Hf 
that is of the order of the typical effective transverse field 
produced by the TD interactions in the system. Note 
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FIG. 3: Spin glass - paramagnet phase transition line in the 
transverse coupling - temperature (A , T) plane, as calculated 
for the Hamiltonian in Eq.lQ (dashed line) and including the 
effective field produced by the off-diagonal dipolar interac- 
tions (solid line). Circles denote experimental data taken from 
Fig.f of Ref. Q. 

that our analysis suggests that it is not sufficient to treat 
the TD interactions in MF, i.e., to replace, e.g., JfJf 
by Jf(Jf) + (J?)Jj. Only the second term potentially 
adds to the effective transverse field, but at the transi- 
tion there is no such effect in this approximation since 
(J?) = 0. 

Finally, we consider the reduction of the cusp in the 
nonlinear susceptibility Xi{T) at low T, found experi- 
mentally by Wu et al. (see Fig. 3 of Ref. 0). This result 
was interpreted as evidence of a first order phase transi- 
tion. However, the experimental results should be recon- 
sidered in view of the significance of the hf interactions, 
leading to Eq.|0, and of the TD interaction. Within the 
model © the transition is expected to be second order. 
However, by taking into account the electro-nuclear na- 
ture of the fsing states and the renormalization of the 
spin, the diminishing of the cusp is naturally explained: 
X3(T) is conjugate to the Edwards- Anderson parameter, 
and is proportional [3 to Y.iA s i. S ]) 2 = 5 oE,,( s ' s j> 2 . 



where So is the effective electronic spin of the state | +) . 
Experimentally, we assume that the transition point can 
be approached with a certain accuracy, giving a maxi- 
mum value to J2ij( s i s j) 2 'i this would result in a maxi- 
mum value Max(x3(T)) cx Sq. As T is reduced H C {T) 
increases, thereby reducing So and Max{xz{T)), until 
the cusp can not be experimentally resolved from the 
background. In order to establish a different scenario 
for the phase transition one must treat properly the TD 
interactions. In Ref. 17] it was shown that a random lon- 
gitudinal field smears the SG-PM transition of the quan- 
tum Ising model, and we expect the TD interaction to 
do the same. 

In this letter we have shown that the usual description 
of the LiHo x Yi_ x F4 system, in which both hyperfine and 
TD interactions are neglected, is not adequate to explain 
the SG-PM phase transition. The large longitudinal hf 
interactions of the Ho ion force the relevant Ising doublet 
to be an electro-nuclear state, with definite and opposite 
values of the spin-7/2 nuclear spins. If it were not for 
transverse hf interactions, the system would then actually 
behave like a classical Ising model with renormalized pa- 
rameters. The quantum nature of the system is revealed 
only at large transverse fields, proportional to the energy 
gap to the first excited electronic state, where quantum 
fluctuations induced by the transverse hf interactions be- 
come significant. The TD interaction serves to enhance 
the effective transverse field, therefore reducing the crit- 
ical field of the transition. A further study of the role of 
the TD interactions would be of interest, since our ap- 
proximate MF treatment does not incorporate multi-spin 
fluctuations Recently many interesting phenomena, 
some connected to the nuclear spins in the system, were 
observed in the LiHo y Yi_^FA[l^ l20ll2I| . Our approach 
to the hf interactions applies to all values of dilution at 
low T , and we hope it will help in understanding these 
phenomena as well. 
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